Optical coupling to spin waves in the cycloidal multiferroic BiFeC>3 
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The magnon and optical phonon spectrum of an incommensurate multiferroic such as BiFeOa 
is considered in the framework of a phenomenological Landau theory. The resulting spin wave 
spectrum is quite distinct from commensurate substances due to soft mode anisotropy and magnon 
zone folding. The former allows electrical control of spin wave propagation via reorientation of the 
spontaneous ferroelectric moment. The latter gives rise to multiple magneto-dielectric resonances 
due to the coupling of optical phonons at zero wavevector to magnons at integer multiples of the 
cycloid wavevector. These results show that the optical response of a multiferroic reveals much more 
about its magnetic excitations than previously anticipated on the basis of simpler models. 



The coupling between different types of correlated elec- 
tron order in "multiferroic" materials possessing simulta- 
neous magnetism and ferroelectricity gives rise to several 
interesting phenomena [H, H & 0] ■ This coupling leads 
to rather complex broken-symmetry phases [5[ and ex- 
citation spectra. Recent examples are the magnetically 
induced ferroelectricity in TbMnC>3 and the switching 
of the magnetic state by an electrical probe observed in 
BiFeOs [6] . Crystals of these two materials show incom- 
mensurate magnetic order: in TbMnC>3, the incommen- 
surate order established at 41 K is believed to generate 
the spontaneous polarization observed at 28 K, while in 
BiFeOs, the cycloidal antiferromagnetic order is estab- 
lished at 650 K, well below the ferroelectric transition 
at 1120 K. The cycloid disappears in a strong (18 T) 
magnetic field [3] or in thin films 0]. 

This paper studies the excitations and electromagnetic 
response in a cycloidal multiferroic such as bulk BiFc03. 
The combination of incommensurate order and magne- 
toelectric coupling is found to give a strong coupling be- 
tween a long-wavelength electric field and spin waves at 
multiples of the cycloid wavevector. This leads to a series 
of resonances in the dielectric constant at integer and ir- 
rational multiples of a fundamental frequency determined 
by the magnon energy at the cycloid wavevector. The 
lowest-lying resonances created by this effect should be 
visible in standard transmission and reflectivity measure- 
ments at room temperature in the far infrared frequency 
range (below the lowest optical phonon at 2 THz). Re- 
cent Raman and optical reflectivity spectroscopy studies 
of bulk BiFeC>3 Q focused on the optical phonon reso- 
nances, without any interpretation of the far IR region. 
The optical phonon spectra seem to be well understood 
from first principles calculations 11]. Nevertheless, the 
latter completely ignores the underlying cycloidal mag- 
netic structure and its magnetoelectric character. Our 
work reveals a rich sub-phonon dielectric response not 
yet explored by optical experiments or first principle cal- 
culations. Aside from allowing a direct observation of the 
basic physics of multiferroic coupling in BiFeOs and re- 
lated materials, this effect could be used in devices based 



on electronic excitation of spin waves [13], or the de- 
velopment of fast magnetic probes for magnetic domain 
switching. 

It has been known for many years that magneto- 
electric coupling in uniformly ordered materials mixes 
spin waves (magnons) and polarization waves (optical 
phonons) [12 . 1 1 31 ] - This gives rise to low- frequency 
magneto-optical resonances in the dielectric susceptibil- 
ity, the so-called electromagnon excitations, that may be 
visible in infrared experiments [1, 0|- However, the reso- 
nance frequencies are primarily determined by the zero- 
wavevector magnon frequencies. The mixing with finite- 
wavevector magnons found here results specifically from 
the incommensurate magnetic order; while such incom- 
mensurate order complicates the calculation of excita- 
tion spectra because a finite-wavevector analysis is now 
necessary, incommensurate order is a common feature 
of many of the most studied multiferroic materials. Al- 
though recent theoretical work has clarified the relation- 
ship between ferroelectricity and spiral ferromagnetic or- 
der 14j, |l5j , the characteristic features of electromagnon 
excitations in incommensurate magnets seem not to have 
been studied before. 

This paper assumes that the ground state consists of 
a uniform polarization and an incommensurate magnetic 
structure, as in many important multiferroic materials; 
an interesting effect when the polarization and magne- 
tization are both periodic and commensurate with each 
other has recently been discussed [lil ]. 

Our calculation is based on a dynamical Ginzburg- 
Landau treatment of the coupled ferromagnetic, anti- 
ferromagnetic, and polarization orders, in a cycloidal 
ground state and an applied AC electrical field. It is 
assumed that the system is far below its critical tem- 
perature so that thermal fluctuations can be ignored in 
the dynamics. The long-period incommensurate cycloid 
is produced in our model by a Lifshitz term in the ef- 
fective Ginzburg-Landau free energy; the existence of 
such a term was previously argued [7, 17, 3, 3 as the 
unique symmetry-allowed explanation for the observed 
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FIG. 1: (color online). Antiferromagnetic lattice with har- 
monic cycloid order. The cycloid plane is pinned by the po- 
larization Po, with arbitrary propagation direction q in the 
plane perpendicular to Po- We also depict the tangent vector 
D related to phase fluctuations of the cycloid. 



order. The model free energy is 
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Here L = \M\ — M 2 \ is a Neel vector describing the 
staggered sublattice magnetization, M = \M\ + M 2 \ is 
the total magnetization of the material, and P z is the 
magnitude of the ferroelectric polarization along z [the 
cubic (111) and equivalent directions in BiFeOa]. Many 
possible terms have been omitted from the free energy 
as absent or unimportant in BiFe03. Note that in the 
absence of the Lifshitz term, the AFM and FE orders 
are decoupled, and the ground state is a commensurate, 
isotropic Heisenberg antiferromagnet with |Lo| 2 = 
and an easy-axis ferroelectric with uniform polarization 
given by P = P z, P 2 = —a/u. Clearly A < and r > 
for an antiferromagnetic ground state, while a < and 
a± > for a ferroelectric ground state. 

The Lifshitz term induces an incommensurate cycloidal 
order in which the antiferromagnetic moment rotates in 
an arbitrary plane including 2. It also increases the mag- 
nitude of Lq and Po- Our reference ground state is a 
cycloid with 



Lq(x) = Lq [cos(gx)z + sin(qa;)x] . 



(2) 



Here the cycloid direction has been chosen along x. The 
pitch is q = aPo/c and the magnitudes of antiferro- 
magnetic and polar order are Lq = (—A + a 2 P^ /c)/G, 
Pq = (—a + o?L\j 'c)/u. A check on the above is that 
there are two zero-energy symmetry actions for fixed P: 
changing the cycloid phase and rotating the cycloid di- 
rection in the xy plane. The displacement for the phase 
change from the transformation qx — * qx + <p/ Lq is 



SL = 4> [cos(qx)x — s\n(qx)i\ 



(3) 



We can use the shorthand D(x) for this spin vector 
direction transverse to the local ground-state spin di- 
rection, and write 5L = <fiT)(x) for the phase change. 
The displacement for the cycloid direction change due 
to an infinitesimal rotation qx — > q(x + r/y) is <5L = 
■qLQqyt){x) + r]L sin(ga;)y. This cycloid direction sym- 
metry is similar to the one found in smectic liquid crys- 
tals in that it requires large displacements relative to the 
original state at some point in the crystal. 

The linearized equations of motion for the ferroelectric 
and antiferromagnetic order parameters are 
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Here 7 is a gyromagnetic ratio with dimensions of (sG) _1 , 
while / has dimensions of s -2 and plays a similar role in 
the ferroelectric equation of motion. In the above po- 
larization equation, we have ignored damping and the 
possible Poisson-bracket term between P and M [20| . 
The coupled spin and polarization equations are solved 
in terms of the parametrization 
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The field cj>(r) denotes phase fluctuations of the cycloid 
ground state, while ip(r) refers to spin fluctuations out of 
the cycloid (xz) plane. The field 5p(r) — P — Pq denotes 
optical phonon fluctuations related to longitudinal and 
transverse vibrations of the ferroelectric moment P [23J . 



Calculating SF/6L to linear order in the displacements, and inserting into Eq. (|4"a|) leads to 
[lj 12 + cV 2 ] <p = 2cqsm(qx)(dyip) - aL (V x Sp) ■ y, 
[a/ 2 - cq 2 + cV 2 ] ip = -2cqsm(qx){d y 4>) - aL 2qcos(qx)(5p y ) - (V x Sp) ■ D 



(6a) 
(6b) 



with u>' — lu/^Lq^/t) defining the frequency in units of ~ 10 12 rad/s. The equation of motion for the polarization 
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becomes 
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where we define £ = r(^L ) \ j f (~ 10 in BiFeOa), aii = —2a + 3a Ig/c, and the AC electric field is i£e 



Some intuition for Eqs. (|6a|) - (I6c|) may be gained by 
considering their invariance under rotations of the cycloid 
in the xy plane. Consider the u) = symmetry operation 
4> ~~ > VQV^ ~ * ?ysin(gx), with <5p left unchanged. This 
creates no terms in both sides of Eq. (fuaj) . In Eq. |6b|) . 
the left hand side becomes — 2cg 2 ?ysin(g;z;), which cancels 
the term generated in the right hand side due to the y- 
dependence in <j>. For the square brackets in Eq. (16c]) , the 
first term adds gcos(ga;)[?7sm(gx)]y, with the third term 
— 'V[rjsin(qx)] x D giving the desired cancellation. 

In the commensurate limit, a, q — > and we recover 
two transverse AFM fluctuation modes with lo' 2 — ck 2 . 
There are two dispersionless polarization wave modes 
transverse to Pq with frequency ui' 2 = a±/£, and one 
longitudinal mode with frequency lo' 2 = ay/£. The com- 
bination orthogonal to <frT), and ipy, S — Lo • SL, is con- 
stant in time in the equation of motion above but decays 
dissipatively if non-Poisson bracket terms are kept. It 
is not expected to modify the above propagating modes. 
For k — ► in the <fr mode, we obtain a symmetry (the 
phase shift) which should be nondissipative. For k — > 
in the ip mode, there should be some finite dissipation. 

We solve the system Eqs. (l6al) - ([6c|) using the ansatz 
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with n an integer running from — oo to oo and the restric- 
tion \k x \ < q/2. The complex numbers (j) n , ip n , {pi) n are 
independent of r but have a k dependence. Substitution 
of this ansatz into Eqs. ([6a|) - ((6c|) yields a linear system 
of equations whose characteristic polynomial determines 
the resonance frequencies. 

Consider the unmixed spin waves in the limit (p) n — > 0. 
Substitution of the ansatz into Eqs. (|6af . (|6bp gives 
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with k 2 = (k x + nq) 2 + k 2 + k 2 z . We see that modes 
propagating along the cycloid plane (with k v = 0) are 
simple plane waves with cyclon (</>) and out of plane (tp) 



dispersions given by lo' 2 = ck 2 and 



c(fc + q 2 ) re- 



spectively. As expected from symmetry, the cyclon mode 
remains soft, but the out of plane mode ip acquires a gap 
due to the pinning of the cycloid plane by the ferroelectric 
moment. Note that the ipo gap equals the cyclon energy 
at n = ±1. At k = 0, the AFM resonance modes are sim- 
ply lj' 71 — yfcq\n\ and ui' n = y/cqV n 2 + 1 for the cyclon 




FIG. 2: (color online). Spin wave spectra for a cycloidal mul- 
tiferroic including one optical phonon mode (not shown) . The 
solid lines are for propagation along the non-trivial direction 
perpendicular to the cycloid plane (y in Fig. 1). The dashed 
lines are the spectra for propagation along z, including zone 
folding. The magnon modes at nq are labelled by <j> n (cyclon) 
and ipn (out of plane). 



and out of plan e excitations respectively. For BiFeOs, 
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eg" ~ 1 |l8l . , and the modes are equally spaced by 
approximately 10 12 rad/s intervals (0.16 THz or 10 K). 

If the propagation vector has a small projection along 
the k y direction, the dispersion curves in the cycloid 
plane will repel each other whenever they intersect. This 
leads to a series of small gaps (anticrossings) in the prop- 
agation frequency as a function of the reduced k x , sim- 
ilar to the effect discussed by Bar'yakhtar in an heli- 
cal ferromagnet [2l|. Apart from these small gaps, the 
high-frequency dispersion curves are nearly unaffected by 
the incommensurate order. A very different situation is 
found for the low-frequency modes along the non-trivial 
y direction. Consider for example the lowest frequency 
mode <fro- Solving Eqs. f[8aj). (|8b|) within 0(k 6 ) leads to 
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This soft mode dispersion is strongly anisotropic, and 
may be useful for electrical control of the spin wave group 
velocity via switching of the Pq direction. A similar effect 
is found for phase fluctuations in smectic liquid crystals, 
and in an helical ferromagnet such as MnSi [22]. The 
anisotropy is related to the q rotation symmetry dis- 
cussed above, which forbids a k 2 term in the dispersion. 

A full numerical solution for propagation along y is 
shown in Fig. [21 Note that in addition to the soft mode 
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anisotropy effect, k y ^ admixes the modes at nq and 
— nq, splitting their dispersions. Moreover, the second 
lowest frequency mode (connected to <j)±x) acquires a neg- 
ative group velocity. 

The full solution of Eqs. (pa|) - ([Sc|) with (Sp) ^ shows 
avoided crossings between the optical phonon dispersion 
and magnon branches at finite k. At these anticross- 
ings the spin wave excitations are highly mixed with the 
polarization wave. However, an optical experiment only 
probes modes at k w due to the large value of the 
velocity of light when compared to the magnon velocity 
(~ 10 6 cm/s). One can ask whether any n > mode re- 
sponds to optical excitation with wavelength much larger 
than 2ir/q ss 600 A. The answer is affirmative as now 
shown by calculation of the AC electric susceptibility. 



An AC electric excitation Ee 



may excite spin 



waves due to the appearence of magnetoelectric fields 
such as 2gcos (qx)(Sp y ) on the right hand side of Eq. j6b|) . 
This couples ip to an optical phonon at k = 0. Solving the 
system of equations shows that the AC electrical suscepti- 
bility tensor (pi)o = XijEj is diagonal for our model, with 
Xxx = -T/(W 2 - a±), Xzz = -1/(W 2 - 0||). Therefore 
the optical response with light polarized within the cy- 
cloid plane contains no magnetic anomaly but only simple 
poles at the bare optical phonon frequencies. However, 
if the light is polarized along the direction perpendicular 
to the cycloid plane, the situation is quite different, 



-(J 2 -2cq 2 ) 
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(10) 



This AC susceptibility has two poles at the follow- 
ing shifted phonon and AFM (ip±i mode) frequencies: 
= OxA + 2cq 2 , ft'l FM = 2cq 2 - 2(aqL ) 2 /a ± . 

Hence we conclude that the simplest model for BiFe03 
already contains a low-frequency magneto-optical reso- 
nance. Note the prominent importance of magnon zone 
folding: Even though the AC electric field is a zero 
wavevector excitation, it couples to a magnon at wavevec- 
tor ±q. This is a special feature of electromagnetic exci- 
tations in an incommensurate multiferroic material. 

The magnon excitations of a cycloid with wavevec- 
tor q may violate momentum conservation by a multi- 
ple of q (a reciprocal lattice vector for the cycloid). In 
the case of a multiferroic material with incommensurate 
magnetic order, this allows the possibility of observing 
additional magneto-dielectric resonances due to the cou- 
pling of a zero-wavevector phonon to magnon modes at 
other multiples of q. To see this effect, consider a contri- 
bution to the free energy in the form of an easy-plane 
anisotropy ^(P-L) 2 , which is known to add anhar- 



monicity to the cycloid order [17(. To first order in A, this 
term adds a third harmonic to Eq. ([2]) . This gives an ad- 
ditional magnetoelectric contribution to the right side of 
the equation determining the cyclon field (p [Eq. (l6a|) ]: 
XPqLq/2 [cos (2qx)(6p x ) — sin (2qx)(Sp z )]. As a result, 
the AC electric susceptibilities Xxx , Xzz develop poles at 



the magnon modes <fi±2- Similarly, the out of plane field 
ip is excited by a term of the form cos (3qx)(Sp y ), show- 
ing that Xyy has an additional electromagnon pole at the 
mode ^±3- Considering additional powers of A (higher 
harmonics of the cycloid) shows that higher frequency 
resonances at nq are present, with strength falling off as 
higher powers of A. 

In conclusion, we determined the electromagnon spec- 
tra for a class of cycloidal multiferroics including bulk 
BiFe03. A simple Landau-Ginzburg model shows that 
a low-frequency optical probe will couple to magnons at 
multiples of the cycloid wave vector, leading to a series of 
magneto-dielectric resonances at integer n and irrational 



multiples of a fundamental cyclon frequency 



and tpn magnons respectively). Interestingly, only one 
mode (ip±i ) becomes electric dipole active due to the 
strong linear magnetoelectric effect driving the harmonic 
cycloid order. The remaining modes acquire a small elec- 
tric dipole character due to weaker quadratic magneto- 
electric effects and anharmonicity in the cycloid order. 
As a result, optical experiments such as transmittivity, 
reflectivity, and Raman spectroscopy will reveal much 
more about the magnetic order of multiferroics than pre- 
viously expected on the basis of homogeneous models. 
For example, the direction dependence and intensity of 
higher resonances is directly related to anharmonicity in 
the magnetic order. Other incommensurate multiferroics 
may show similar phenomena even if the details of the 
multiferroic coupling are different. 
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